Fractional vortex in asymmetric 0-n long Josephson junctions 
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We consider an infinitely long 0-7r Josephson junction consisting of and tt regions having dif- 
ferent critical current densities j Ct o and j Ctn . The ground state of such a junction corresponds to 
a spontaneosly formed asymmetric semifluxon with tails decaying on different length scales. We 
calculate the depinning current of such a fractional vortex and show that it is different for positive 
and negative bias polarity. We also show that upon application of a bias current, the fractional 
flux (topological charge) associated with the vortex changes. We calculate the range of fractional 
flux associated with the vortex when the bias changes from negative to positive critical (depinning) 
values. 



in 

(N 



PACS numbers: 74.50.+r, 85.25.Cp 

Keywords: 0-tt long Josephson junction, fractional vortex, semifluxon 



C3 
O 
O 

Oh- 
33 



-a 

C3 

o 
o 



> 
o 
00 
o 

rn 
O 
m 



x 



I. INTRODUCTION 

A 0-tt Josephson junction (JJ) combines the properties 
of a conventional J J with critical current density j c _o > 
with the properties of a tt J J with j c> „ < 0. Such JJs 
play an important role in the determination of the order 
parameter symmetry in novel superconductors^— , allow 
experiments with fractional Josephson vortex matter— ~— 
and recently were used as a way to construct a tp JJ with 
a tunable current-phase relation^. 

Nowadays there are several technologies that allow to 
fabricate 0-tt long Josephson junctions (LJJs)^^^— . 
One of themiSr— is based on employing a ferromagnetic 
barrier, which has a different thickness of the ferromag- 
net in the and tt part and, correspondingly, different 
critical current densities j c .o and j c ,w I n practice, it is 
very difficult to control the thicknesses of the barrier in 
the and tt region very precisely. Therefore, j' c and j Ci7r 
are always different by absolute value, i.e., j c ,o 7^ |j c .ir|- 
Moreover, for some devices such as tp JJs^ or, more gen- 
eral, for JJs with a tunable current-phase relation&ii, 
it is even necessary to make j Cj o and |j c ,ir| different to 
achieve the required properties. However, most of the 
theoretical works so far deal with the idealized situation 

jc,o = \ic,A- 

Therefore, in this paper we consider an infinitely long 
0-7r J J with jc.o 7^ |jc.7r| and investigate the qualitative 
differences in comparison with the symmetric case. 

The paper is organized as follows. In Sec. [II] we inves- 
tigate the ground state and show that it corresponds to a 
semifluxon (fractional Josephson vortex carrying the flux 
±$o/2), which is pinned at the 0-7r boundar y 15 : 16 , but 
has asymmetric tails. At non-zero bias current the shape 
of the fractional vortex can be obtained only numerically. 
Similar to the case of a symmetric 0-tt LJJ, by applying 
a large enough bias current, one can reach the depinning 
current^! - — , investigated in detail in Sec. IIIII We will 
see that in an asymmetric 0-7T JJ the critical current is 
different for positive and negative bias polarity. More- 
over, we identify two different mechanisms of switching 



to the non-zero voltage state: depinning of the fractional 
vortex and primitive switching, i.e., when the bias cur- 
rent density exceeds j Ct o or |i c ,7r|- Further, in Sec. IIVI we 
show that the fractional magnetic flux localized at the 0- 
tt boundary can deviate from ±$o/2 when a bias current 
is applied. Finally, Sec. IVl summarizes this work. 



II. FRACTIONAL VORTEX IN AN 
ASYMMETRIC O-tt LJJ 

Consider an infinite asymmetric 0-tt LJJ with criti- 
cal current density j Ci o > in the half (x < 0) and 
jc.Tr < in the tt half (x > 0). For theoretical analy- 
sis it is convenient to work in normalized units. There- 
fore we introduce normalized critical current densities 

70 = ic,o/jAbsAv and 7^ = J CjW /iAbsAv with jAbsAv = 

(ic,o + |jc.7r|)/2 > 0. Then the normalized critical current 
current along the 0-tt LJJ is given by 



3c(x) 



70 > x < 
7^ < x > 



(1) 



A normalized Ferrcl-Prangc equation describing the 
static solutions for the Josephson phase <fr(x) reads 



4>xx - j c {%) sine/) = -7. 



(2) 



Here 7 is the bias current density normalized to jAbsAvj 
subscripts xx denote the second derivative with respect 
to x, while the coordinate x is normalized to 



Aj(ja 



bsAvJ 



27r^ d'jAbsAv ' 



(3) 



calculated using the critical current density jAbsAv In 
Eq. ([3]) the quantity hqcL' is the inductance (per square) 
of the superconducting electrodes forming the LJJ. 

For further analysis it sometimes will be convenient to 
introduce the deviation 8 of 70 and 7^ from the symmet- 
ric case (±1 in our normalized units), i.e., 

la = 1 + 8, i w = -l + 8. (4) 
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Thus one can use a single asymmetry parameter \6\ < 1 
instead of 70 and j w . 

For 7 = the semifluxon solution of Eq. ([2]) is given 
by the fluxon tails properly jointed at the 0-tt boundary, 



III. CRITICAL CURRENT IN AN 
ASYMMETRIC O-tt LJJ 

A. Depinning vs. primitive switching 
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x < 



x > 



(5) 



where An = l/-y/7o and A^ = l/vh^l are normalized 
local Josephson penetration depths in the and tt parts. 
They are given in units of Aj(jAbsAv)- The phase 0(0) 
and its derivative (oc magnetic field) <fi x (0) must be con- 
tinuous at x = 0, i.e., 
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where 



Z§ = cxp 



Ao 



cxp 



Ayr 



By taking tan(. . .) of both sides of Eq. ^ we get 

z n - 1 



zo = 



z-k + 1 



(6) 
(7) 

(8) 
(9) 



Solving Eqs. ([7]) and © for z and z^ we obtain two roots 
(for each of them). The negative roots can be neglected 
since zq and z v are positive by definition ([5]). Thus, the 
remaining positive roots are 




> 1; 



< 1. 



(10a) 
(10b) 



From here, taking into account definitions ([5]), we get 



xo 




Hll < 0; (11a) 
7o / 



7o 



>0. (lib) 



For 7/O the static solutions of Eq. @ cannot be 
obtained analytically. Therefore, Fig. [T] shows profiles 
of the phase 0(x), magnetic field 4> x (x) and supercurrent 
j c {x) sin 4>{x) obtained numerically. For zero bias current 
the profile coincides with the one given by Eq. ([5]) with 
Xq and x v given by Eq. (fTTj) . 



Althouh for 7 ^ the analytical solutions for the frac- 
tional vortex shape cannot be obtained, it is still possible 
to obtain the range of 7 where such a static solution ex- 
ists. By multiplying Eq. ([2]) by 2cf> x and integrating, one 
finds (separately for and 7r region) 



± I V 2 Wo - 7o cos 4> - 70] x<0 ^ 
\ ^f2 [C' n - -f n cos (j) - 70] x > ' 



where the integration constants Co and C w can be found 
from the boundary conditions at 



(13) 
(14) 
(15) 



Conditions (|T4|) and (|T5|) assume a vortex of positive po- 
larity. For Co and C w one gets 



: (±oo) = 0; 

i(-oo) = arcsin(7/7o); 

i(+oo) = 7r — arcsin(7/7 ?r ). 



y arcsin ( ^ ) + V^o ~ 7 



(7) 



tt — arcsin 



(16) 

TtI - 7*. (17) 



In Fig. [2] we plot the phase plane curves 0^(0) given by 
(|T2l separately for the and the 7r region and for different 
values of the bias current 7. At 7 = 0, see Fig. [2ja), the 
semifluxon solution corresponds to a line on the phase 
plane 0-0^ starting at = arcsin(7/7o), <p x = at x = 
—00 (black dot), and going towards the crossing point 
of black and gray curves. At the crossing (half-black- 
half-gray dot) at x = we switch to the gray 7r-region 
curve and follow it up to the point = tt — arcsin(7/7 7r ), 
4>x = (gray dot). Upon increasing 7 the crossing point 
of the two trajectories above transforms into a touching 
point at the critical bias current density j c+ > 0, see 
Fig- Hlb). After this the trajectories disconnect and the 
static solution is lost. Similar things happen for negative 
bias current density, see Fig. [^c) shown for 7 C _ < 0. At 
the touching point, i.e., at = 0(0) and at j c (it can be 
either 7 C+ or 7 C _), the <f> x as well as d(4> x ) / d<f) of both 
trajectories are equal, i.e., 

Co (7c) ~ 70 cos 0(0) - 7c0(O) 
= C T ( 7c )-7 w cos0(O)-7 c 0(O); (18) 
70 sin 0(0) = 77r sin 0(0). (19) 

From Eq. (fT9"l) we conclude that either 0(0) = tt (for 
7 C > 0) or 0(0) = (for 7c < 0). Then Eq. (HHJ) becomes 
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FIG. 1: (Color online) The phase 4>{x) (a), the magnetic field <f> x (x) (b) and the supercurrent jo,* sin<j)(x) (c) of a postive vortex 
in a 0-7r LJJ with the critical current density asymmetry 8 = 0.2 (left pannel), corresponding to depinning, and 8 = 0.4 (right 
panel) , corresponding to promitive switching. Profiles are shown for zero bias, positive precritical bias, and negative precritical 
bias. The gray areas in (c) indicate the critical current densities ±70 and ±7^ in and it parts. 



n (70) + arcsin (t 1 ) + ~ ^ ~ V ^ _= 7c2 = nic T ~ 



(20) 



r 



where the upper sign corresponds to the positive j c = 
7c+ > and the lower sign to a negative 7 C = j c - < 0. 

Note that if 70 = — 7^, so that the l.h.s. of Eq. (|2TJ)) 
vanishes, one obtains the well known depinning current 
of a semifluxon 7 C = ±2/n. For the asymmetric case 
70 7^ —7^ Eq. (|20[) provides an implicit dependence of 



7 C on 70 and 7^. Note, however, that for given 70 and 
7,r, 7 C must be sought in the domain |7c| < min(7o, |7tt|), 
otherwise some of the square roots and arcsin functions 
in Eq. (|20p are undefined in the real domain. 

In terms of 5 Eq. (|20p reads 
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+ y/(l + 8f - 7 2 - y/(l - (5)2 - ^ = ^ _ 2 Sgn( 7c ), 



(21) 




equation for <5 C 
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and 7 C should be searched in the domain | 7c | < min(l + 
5,1 — 6). One may or may not have solutions in this in- 
terval. And the solution of Eq. (f2Tj) (dis) appears exactly 
when 



\-y c {6)\ = min(l + 6, 1 - 6). 



(22) 



Thus, we can find a "critical" value 6 C when solution(s) 
(dis)appear. 

First, for 7 C > and assuming 6 > we substite j c = 
1 — 6 from Eq. ([2"2]) into Eq. (|2"Tj) and get the following 
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(1 - 6 C ) + 2y/6 c + 2 = 0, (23) 



which can be solved numerically to give 6 C w 0.2606. For 
6 < 0, we substitute j c = 1+6 from Eq. (JUJ) into Eq. (|2lj) 
and obtain an equation, which has only the trivial solu- 
tion 5 = — 1. Thus, for 7 C > Eq. ([2"Tj) has a single 
solution only in the interval 5 = — 1 . . . 6 C . 

Second, for 7 C < 0, similarly to the previous case we 
obtain that Eq. (|2ip has a single solution, if S = —5 C ...+ 
1. 

Finally, when the asymmetry 8 is such that Eq. (|2Tj) 
delivers no solution for 7 C (no vortex depinning current) 
for positive or negative bias polarity, the switching to 
the voltage state takes place at 7 C given by Eq. (|22|) . i.e., 
when the critical current density in the or 7r part will 
be exceeded. This can be written as 



7c+ 



V- 

and, similarily. 



fPositiveSolutionOfdH]) for - 1 < 5 < 8 C 



for 6 C < 5 < +1 



(24) 



7c- = 



-1-6 

NegativeSolutionOf 



for - 1 < 5 < —6 C 
for - 6 C < 6 < +1 ' 
(25) 

Fig.[3]shows the resulting dependence 7 C (£) for positive 
and negative currents. One can see that in the general 
case (5^0 the positive and negative critical (depinning) 
currents are not equal by absolute value. For 8 = ±8 C 
the pinning has a maxumim in one of the bias directions. 
For 1 6 1 = 1, i.e., when one of the two halves has no 
Joscphson properties (zero critical current density) the 
critical current vanishes. 

It turns out that one can find a very good analytical 
approximation to the solution 7 C (<5) of the transcendental 
Eq. (JUJ. If 8 < 1, one can solve Eq. j3T]) by Taylor 
expanding it near 8 = up to terms ~ 8. One obtains 



7r 7c T 2 = 2Sy/l - 7 2 

The solution of this equation is 



7c± 



sVW 2 



TT 2 ± 7T 



4<5 2 



(26) 



(27) 



The formula (|27p gives a very good approximation to 
the 7c (<5) dependence obtained numerically by solving 
Eq. (|2"Tj) . as can be seen in Fig. [3] 

Finally, we have checked our results by numerically 
solving a timc-dependent sine-Gordon equation using 
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FIG. 3: (Color online) The dependence of the normalized 
critical current 7 C ± of an asymmetric 0-n LJJ on the asym- 
metry parameter S. The fat black continuous line shows the 
final 7 C ± (S) dependence for all 5 given by Eqs. (|24[) and (|25|l . 
It consists of vortex depinning branches (curved) given by 
Eq. (|21[) and primitive switching branches (thin dash, straight 
part of fat black continuous line) given by Eq. (|22p . The ap- 
proximation Eq. (|27p is shown by the fat dashed line. Gray 
curves marked as — r y c -{8) and — j c +(8) represent the critical 
current of the fractional antivortex. 



StkJJ2£. A perfect agreement between the theory de- 
velopped above and direct numerical simulation is found, 
including a crossover between different branches. Numer- 
ical results are not shown in Fig. [3] to avoid overcrowd- 
ing. In Fig.[T]we also show the profiles of the phase (f){x), 
phase gradient 4> x (x) (oc magnetic field) and supercurrent 
j c (x) sin (f>(x) at precritical bias currents density. One 
can see that just before depinning (6 = 0.2 and j c > 0) 
the second dcriative of the phase (supercurrent) becomes 
continuous and has the same sign everywhere. On the 
other hand, just before primitive switching (5 = 0.4 and 
7 C > 0) the discontinuity of the supercurrent persists. 
These precritical profiles of the supercurrent are exactly 
the ones that should be seen in supercurrent distribution 
images obtained by low temperature scanning electron 
microscopj^^i. 

The numerical simulations also reveal the branch cor- 
responding to the critical current of a fractional antivor- 
tex. Indeed, the theoretical description presented above 
applies to a fractional vortex with positive flux (topolog- 
ical charge), sec Eqs. (|T4l) and ([151) . For an "antivortex" 
with negative flux (polarity) the situation reduces to the 
one considered above if one inverts the sign of the bias 
current. Thus, for an antivortex, the positive critical 
current is —j c -(6), while the negative one is — 7 c +(5). 
These curves are shown in Fig. [3] as well. One can see 
that, e.g., for 6 < 0, the positive vortex becomes un- 
stable at 7 C _|_((5). It emits an integer fluxon and turns 



into a fractional "antivortex" , which is still stable at this 
value of bias current. Depending on parameters such as 
damping, one will observe either the critical current of 
a fractional vortex 7 C+ (<5) or the one of the "antivortex" 
-7c- (S). 



B. Dynamics of depinning 

We have studied numerically the dynamics of the 
switching to the resistive state and compared the case 
of depinning at 6 = 0.2 < 6 C with the case of primitive 
switching at 6 = 0.4 > S c . 

For 6 < 5 C the depinning at 7 C + starts from flipping 
the fractional vortex into an "antivortex" and emission 
of a fluxon. Further dynamics depends on the stability 
of the "antivortex" at given bias. First, if the "antivor- 
tex" is unstable, which is the case if its critical current 
— 7 C _ < 7c+, i.e., at 5 > 0, then the flipping of the frac- 
tional (anti)vortex and the emission of fluxons and anti- 
fluxons continues. This type of depinning dynamics was 
discussed earlier for symmetric LJJsil. Second, if the 
"antivortex" is stable, which is the case for 5 < 0, then 
further flipping does not take place and the state of the 
system is defined by the destiny of the emitted fluxon. 
For large damping the emitted fluxon moves away and 
is absorbed at the far end of the LJJ so that (after this 
transient) the LJJ remains in the static situation with 
an "antivortex" trapped and one finally measures in ex- 
periment/simulation the critical current — 7 C _(<5) of an 
antivortex. For moderate damping the emitted fluxon 
reflects from the edge of the LJJ as an antifluxon and 
starts moving back and forth (colliding with the frac- 
tional "antivortex" near x = on every pass) similarly 
to the dynamics at a zero-field step so that one detects 
a critical current j c +(,6) and a step on the IVC at larger 
bias current, presumably up to —^ c -{6). 

For 5 > S c , depinning at 7 C + starts from switching the 
"weaker" (^l < 70) 7r-part of the LJJ into the resistive 
state, which may look like the penetration of an avalanche 
of fluxons from the edge of the LJJ. One switches to a 
finite voltage state and detects j c +(S) as a critical cur- 
rent. Similar types of switching dynamics is observed at 
the negative critical current. 



IV. FLUX LOCALIZED AT A O-tt BOUNDARY 

From Eq. (|14p and fT5|) it follows that the topological 
charge of the vortex is equal to ir (by absolute value) only 
if 7 = or 70 = \"fn\- Otherwise it is given by 



p( 7 ) = 7r - arcsin(7/7 7r ) - arcsin(7/7o) 

7 \ • ( 7 



7T — arcsm 



■1 + 6 



arcsm 



1 + 6 



(28) 



Naively, from Eq. (|28[) . one concludes that in the case 
when |7 W | <C 70 (6 —> 1) one has p — s- 3ir/2 at 7 —> 7^ = 



6 



(5 — 1. In the opposite case |7 ff | ^> 70 (<5 — > — 1) one 
has p— ?>7r/2at7^7o = 1 + 5. However, this is not 
completely true. 
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FIG. 4: (Color online) The topological charge of the vortex 
p at positive and negative critical current -y Ct ± as a func- 
tion of asymmetry parameter <5 calculated numetrically using 
Eq. ([28j) with 7 = 7 C given by Eqs. <[2TJ) and ([22j). 

Since the topological charge given by Eq. (f2"5)l is a 
monotonous function of 7 it reaches its extremum at 
7 = 7 C ±. In Fig. [4] we plot the topological charge p (the 
flux $ = 3>op/27r) reached at the positive and negative 
critical current 7 C ± as a function of asymmetry d. When 
<5 = (symmetric case) the topological charge p = n 
does not change with bias (see the crossing point of two 
curves). If <5 7^ 0, the topological charge p = it at 7 = 
(horizontal dashed line) and changes in the interval be- 
tween two curves upon application of bias. Each of the 
p(7 c= t) curves has a breaking point at ±<5 C , corresponding 
to a crossover from vortex depinning to primitive switch- 
ing, see Fig. [3] One can see that, indeed, the primitive 
switching branch gives the maximum topological charge 
pmax = 37t/2 ($ max = 3$o/4). However, the depinning 
branch never reaches p = tt/2. To calculate the minimum 
topological charge, one has to derive the limiting behav- 
ior of 7 C+ at 6 — > — 1. By Taylor expanding Eq. f|2 1[) near 
<5 = — 1, we find that 

7 C+ a(5+ 1) at S -> -1, (29) 

where a ss 0.7246 is a solution of the trancendental equa- 
tion 

1 + yl — a 2 + (arcsin a — ir)a = 0. (30) 

Substituting the expression (|29| instead of 7 into Eq. (|28|) 
for (5 — > — 1 we obtain 

Pmin = 7T — arcsin(a) ps 0.74187T. (31) 



This corresponds to $ m i n ps 0.3709$o and is shown as 
the horizonal (short dashed) line in Fig. |4j 

At the end we would like to point out that the topo- 
logical charge of an "antivortex" is, similar to Eq. (|28jl. 
given by 

p = — 7r — arcsin (7/7^) — arcsin (7/70) 

= — 7r — arcsin ( — - — - ) — arcsin ( „ ^ . ) .(32) 

One can see that it is opposite to p of the positive vortex 
only if 7 = or 70 = |7^| (<5 = 0). In the general case 
the topological charges (fluxes) of positive and negative 
fractional vortices at a given bias current are not equal. 
That is why we put the word "antivortex" in quotation 
marks throughout this paper. 



V. SUMMARY 

We have demonstrated that an infinitely long 0-7T LJJ 
with unequal critical current densities 70 7^ {j^l in and 
7T parts has a ground state corresponding to an asym- 
metric vortex of supercurrent carrying the magnetic flux 
±$ /2. The tails of this vortex decay on different length 
scales cx 1 / v /C 7o an d oc 1 / ^/[-y^ | as one goes away from the 
0-7T boundary. Upon application of a small uniform bias 
current the vortex deforms but does not move away under 
the action of the Lorenz force as it is pinned at the 0-7T 
boundary. The fractional flux associated with the vortex 
in this state differs from ±<f>o/2, as given by Eq. (|28|) . 
This situation persists up to the critical bias current, 
which is different for positive and negative bias polarity, 
see Eqs. ((24)) and (f2"5]) . Two mechanisms of switching 
to the resistive state are identified: (a) depinning of the 
vortex at the depinning current given by Eq. (f2"Tj) and 
(b) primitive switching when the bias current 7 exceeds 
70 or 7^ given by Eq. (|22|) . The maximum deviation of 
the fractional flux from 3>o/2 is reached at the positive 
and negative critical currents. For the most asymmetric 
0-7T junction (e.g. 7^ = 0) the minimum possible flux 
$min ~ 0.3709$Oj while the maximum possible flux is 
$ max = 3$ /4. Finally, the topological charges (mag- 
netic fluxes localized near the 0-7T boundary) of positive 
and negative fractional vortices at a given bias current 
are not equal, except for the case of zero bias current. 
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